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^ ■ Abstract 



We compute the large scalar four-point correlation functions in general single field inflation 
models, where the infiaton Lagrangian is an arbitrary function of the infiaton and its first 
derivative. We find that the leading order trispectra have four different shapes determined 
by three parameters. We study features in these shapes that can be used to distinguish 
among themselves, and between them and the trispectra of the local form. For the purpose 
of data analyses, we give two simple representative forms for these "equilateral trispectra" . 
We also study the effects on the trispectra if the initial state of inflation deviates from the 
standard Bunch-Davies vacuum. 
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1 Introduction 



Primordial non-Gaussianity is potentially one of the most promising probes of the inflation- 
ary universe [1]. Like the role colliders play in particle physics, measurements of primordial 
non-Gaussian features provide microscopic information on the interactions of the infiatons 
and/or the curvatons. Constraining and detecting primordial non-Gaussianities has become 
one of the major efforts in modern cosmology. Theoretical predictions of non-Gaussianities, 
especially their explicit forms, play an important role in this program. On the one hand, 
they are needed as inputs of data analyses [2-5] which eventually constrain the parameters 
defining the non-Gaussian features; on the other hand, different forms of non-Gaussianities 
are associated with different infiaton or curvaton interactions, and so if detected can help us 
understand the nature of inflation. 

A variety of potentially detectable forms of non-Gaussian features from inflation models 
have been proposed and classified, in terms of their shapes and running. The scalar three- 
point functions, i.e. the scalar bispectra, are by far the most well-studied. For single field 
inflation, a brief summary of the status is as follows. Minimal slow-roll inflation gives 
undetectable amount of primordial non-Gaussianities [6-8]; non-canonical kinetic terms can 
generate large bispectra of the equilateral shapes [9, 10]; non-Bunch-Davies vacuum can 
boost the folded shape [9,11,12]; and features in the Lagrangian (sharp or periodic) can give 
rise to large bispectra with oscillatory running [13,14]. Multifield inflation models provide 
many other possibilities due to various kinds of isocurvature modes, such as curvatons [15], 
turning [16-21] or bifurcating [22, 23] trajectories, thermal effects [24, 25] and etc. These 
models give many additional forms of large bispectra, notably ones with a large local shape. 

We will be getting much more data in the near future from new generations of experi- 
ments, ranging from cosmic microwave background, large scale structure and possibly even 
21-cm hydrogen line. Compared with the current WMAP, these experiments will be mea- 
suring signals from shorter scales and/or in three dimensions. Therefore a significant larger 
number of modes will become available. This makes the study of four- or higher point 
functions interesting, as they provide information on new interaction terms and refined dis- 
tinctions among models. In this paper we extend the work of Ref. [9] and classify the forms 
of large scalar trispectra (i.e. the scalar four-point function) in general single field infla- 
tion models. There have been some preliminary works in this direction [26,27], calculating 
contributions from the contact interaction diagram (Fig. [1] (A)). For models with a large 
trispectrum, there is yet another set of diagrams involving the exchange of a scalar (Fig. [1] 
(B)) that contributes at the same order of magnitudell] In this paper, we complete this 

^Note that for slow-roll inflation, the contribution from the scalar-exchange diagram Fig. [T] (B) is sub- 
leading, while the graviton-exchange contribution belongs to the leading order [28] . 
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(A) 



(B) 



Figure 1: Two diagrams that contribute to the large trispectra. 



program and classify all possible shapes arising in this framework. 

For the bispectra in general single field inflation, the leading large non-Gaussianities have 
two different shapes controlled by two parameters [9]. As we will see here, for trispectra, we 
have four different shapes controlled by three parameters. Some of them have complicated 
momentum-dependence. For the purpose of data analyses, we give simple representative 
shapes that can capture the main features of these functions. We point out the features in 
the shapes that can be used to distinguish among themselves, as well as to distinguish them 
from the trispectra of the local form. We also study the effects of a non-Bunch-Davies initial 
state of inflation on these trispectra. 

This paper is organized as follows. In Section 2, we review the basic formalisms and main 
results for the power spectrum and bispectra in general single field infiation. In Section 3, 
we calculate the leading order trispectra, and summarize the final results. At leading order, 
the trispectra can be classified into four shapes, controlled by three parameters. In Section 
4, we investigate the shapes of the trispectra, including consistency relations, figures in 
various limits, and also give two simple representative forms of these equilateral trispectra 
to facihtate future data analyses. In Section 5, we discuss DBI and K-inflation as two 
examples to illustrate our results. In Section 6, we study the trispectra when the initial 
state of inflation is in a non-Bunch-Davies vacuum. We conclude in Section 7. 

2 Formalism and review 

In this section, we review the formalisms and main results of Ref. [9]. As in Ref. [29], we 
consider the following general Lagrangian for the inflaton field 0, 




(2.1) 



where X = —-^ 



\g'^'^dn(f)dy(j) and the signature of the metric is (—1,1, 1,1). 
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Irrespective of the specific mechanism that is responsible for the inflation, once it is 
achieved we require the following set of slow-variation parameters to be small, 

where H is Hubble parameter and 

^ P,x + 2XP,^^ ^^-^^ 

is the sound speed. The slow- variation parameters can be large temporarily or quickly 
oscillating [13, 14,30], but we do not consider such cases here. 

The power spectrum is defined from the two-point function of the curvature pertur- 
bation C, 

(C(ki)C(k2)) = (27r)553(ki + k2)^Pc • (2-4) 

For the class of inflation models that we consider, 

1 



(2.5) 



In order to parametrize the three-point function, we need to define a parameter A/E 
related to the third derivative of the inflaton Lagrangian P with respect to X, 

A = X^P,xx + lx^P,xxx , (2.6) 

E = XP,x + 2X^P,xx = ^ . (2.7) 

The bispectrum form factor A{ki, k^-, kz) is defined as 

(C(kOC(k2)C(k3)) = (27r)^53(ki + k2 + k3)P^2np ^ . (2.8) 

i=l '^i 

Up to C(e), this bispectrum is determined by five parameters, c^, A/S, e, r] and s. For the 
most interesting cases <S 1 or A/E ^ 1 where the non-Gaussianities are large, the leading 
bispectrum is given by 
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So we have two different forms determined by two parameters. In such cases, the effect of 
the non-canonical kinetic terms of the infiaton has to become large enough so that the infla- 
tionary mechanism is no longer slow-roll (in slow-roll the canonical kinetic term dominates 
over the non-canonical terms). Since inflation gives approximately scale-invariant spectrum, 
ignoring the mild running of the non-Gaussianity [31], the bispectrum is approximately a 
function of two variables in terms of the momentum ratios, /c2/^i and k^/ki [32]. The two 
forms in (12. 9p have very similar shapes and they are usually referred to as the equilateral 
shapes. Because the two shapes do have a small difference, for fine-tuned model parameters 
Cs and A/S, they can cancel each other to a large extent and leave an approximately orthog- 
onal component. One can use this component and the one of the originals to form a new 
bases of the shapes! 



3 Large trispectra 

As in the bispectrum case, we are most interested in cases where the trispectra are large. 
In general single field inflationary models, this is achieved by the non-canonical kinetic 
terms. The origin of large non-Gaussianities come from terms with derivatives of the infiaton 
Lagrangian P with respect to X . The contribution from the gravity sector is negligibly small. 
The derivative of P with respective to is also small due to the approximate shift symmetry 
associated with the infiaton. Another equivalent way to see this is to work in the comoving 
gauge [6] where the scalar perturbation only appears in the metric. So explicitly does 
not appear in the expansion. Using the leading order relation 

C^-f« (3.1) 

to convert C into a = 50, again we see that P<^ does not appear. 

Therefore for our purpose, it is convenient to choose the infiaton gauge where the scalar 
perturbation only appears in the infiaton [6], 

= 0o(t)+«(t,x) ; (3.2) 

and when we expand the infiaton Lagrangian P, we only concentrate on terms that have 
derivatives with respect to X. Such a method has also been used in Ref. [34,35]. 

^ We would like to thank Eiichiro Komatsu for discussions on this point [33] . 
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3.1 Scalar-exchange diagram 



In this subsection, we compute the scalar-exchange diagram, Fig. [T] (B). Using the inflaton 
gauge, we get the cubic terms of the Lagrangian in the small Cg or large A/S limit, 



2 ' " 6 
Written in terms of C using (13. ip . we get 



A = -2a^^C^ + a|j(l-aC(9.C)^ 



(3.3) 



(3.4) 



Despite of its different appearance from the three leading cubic terms in [9], one can show, 
using the linear equation of motion and integration by part, that the difference is a total 
derivative. 

In terms of the interaction Hamiltonian, = —£3, we denote the two terms in (13. 4p as 



where 



(3.5) 



(3.6) 



Ha{t) = 2a'^Jll^^UPut)UP2,t)UP3,t)i2nr^^^ 

^ 1=1 ^ ^' 1=1 

m) = /n^(p2-P3)6(Pi,t)C/(P2,t)0(P3,t)(27r)35^(X^p.)(.3.7) 

^ i=l ^ i=l 

The is in the interaction picture and satisfies the equation of motion followed from the 
kinematic Hamiltonian. 

The scalar trispectrum is the expectation value of the curvature perturbation in the 
interaction vacuum. According to the in- in formalism [36], there are three terms contributing 
to the diagram Fig. [T] (B)|f 



(C^) = (o| 



ij!dt'Hj{t') 



C/(ki,t)C/(k2,t)C/(k3,t)C/(k4,t) 



if!dt'Hiit') 



|0) (3. 



D fdt' f dt" (0| Hiit') Ci Hi{t") |0) 
ft j-t' 

dt' / dt" (0| Hi{t") Hi{t') Cf |0) 

in ^ to 

dt' / dt" (0| Cj Hi{t') Hi{t") |0) . 

to J to 



(3.9) 



■^The in-in formalism is often used in the literature in terms of a commutator form which is equivalent 
to the form presented here. However, for a subset of terms, the algebra in the commutator form is more 
complicated than the one we use here. We discuss this equivalence in Appendix VK\ 
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Here t is a time several efolds after the modes exit the horizon and to is a time when modes 
are all well within the horizon. In terms of the conformal time r, dt = a{T)dT, we take r = 
and To = — oo. 

We evaluate fl3.9p using the standard technique of normal ordering. We decompose 
(omitting the subscript "I" for ( in the following) 

C(k, t) = C + C = uiK r)ak + w*(-k, T)al^ , (3.10) 

where 

u{k, r) = -7^=il + tkc,T)e-^'^^^ , (3.11) 

and 

[a^,al] = {27if6%k-p) . (3.12) 

After normal ordering, the only terms that are non-vanishing are those with all terms con- 
tracted. A contraction between the two terms, C(k, t') (on the left) and ({p,t") (on the 
right), gives 

[eik,r'),Cip,r")]=uik,r')u*i-p,r"){27r)'5'ik + p) . (3.13) 

We sum over all possible contractions that represent the Feynman diagram Fig. [T](B), where 
the four external legs are connected to C(ki,t)'s. 

To give an example, we look at the 1st term of fl3.9p with the component fl3.6p . One 
example of such contractions is 

c(Pi, toc(p^(i ^3, toc(k'i, t)c(kl, t)c(k3, t)c(k4, t)^i, nai, fxig^ ^, f) 

= [C+(Pi,0,r(ki,t)][C+(P2,0,r(k2,t)][C''(k3,t),r(qi,O][C+(k4,t),r(q2,O] 
[C(P3,0,r(q3,n] ■ (3.14) 

There are three ways of picking two of the three pj's (q^'s), so we have a symmetry factor 
9. Also, there are 24 permutations of the kj's. The overall contribution to the correlation 
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function is 



9-4— ul{t)ul{t)uk,{tW,{t) 
t pt 



dt' / dt"a' {t')a' {t")u,, {t')u,,, {t')ul {t")ul {t")u,,, {t')ul^ {t") 

to J to 



X {27rf6^{Y^ ki) + 23 perm. 



i=l 



9 /A\' k 



12 



1 



8 \Sy kik2kzki {ki + k2 + ku)^ {ks + k4 + k 



12 



{27TfP^6'\J2^i) + 23 perm. , 



(3.15) 



i=l 



where 



^12 



ki + ka 



(3.16) 



The 2nd and 3rd term in (13. 9p has a time-ordered double integration, and so is more 
comphcated. Their integrands are complex conjugate to each other, and we get 



2 ■ 



9 /A 



k 



12 



8 vSy kik2ksk4 {ki + k2 + ku) 



X 



\, iki + k2 + kuf , ^ ki + k2 + ku , 1 



(27r)^P^35^(^ ki) + 23 perm. 



(3.17) 



i=l 



where 



K = ki + k2 + k3 + k4 . (3.18) 
The other terms are similarly computed. We leave the details to Appendix [Bl 



3.2 Contact-interaction diagram 

In this subsection, we compute the contact-interaction diagram. Fig. [T](A). We define 

fi = ]^X^P,xx + 2X^P,xxx + \x^P,xxxx . (3.19) 

^The integrations are conveniently done in terms of the conformal time r. Integrals such as 
j'^^dx x^e^^^ — ±2i are constantly used in the evaluation in this paper. As in [6], the convergence at 
X — cxD is achieved by 2; ^ x{l =F ie). 
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The fourth order expansion is [26] 



Generally speaking, the Lagrangian of the form 

A = ^7oC' + ^7iC' + ^72C + J3 , 
gives the following interaction Hamiltonian at the fourth order in (j [26], 



(3.20) 

(3.21) 
(3.22) 
(3.23) 



4/0 / \ Jo 



9i92 
fo 



9l 



-MO.--. 



(3.24) 



where /, g, h and j's are functions of (, di( and t, and the subscripts denote the orders of 
C. So for (Km we have 



ni = ^(-/. + 4)Cf + |j(3Ac^-s(i-a)(90)^Cl 



(3.25) 



Note that in the second term the order A term cancelled, in the third term the order E term 
cancelled. 

The following are the contributions to the form factor T defined as 

(0 = (2.)^p3^3(^k.) nlr. 



1=1 



i=l 



(3.26) 



The contribution from the first term in fl3.25p is 



36(^-^ 



fc?A:|(k3-k4) 
^3 



from the second term, 

1 /3A 1 ^ 

from the third term, 

1 A (ki.k2)(k3-k4) 

32 Ic? 



\ ^ 3(fc3 + ki) ^ I2hh 



K 



(3.27) 



+ 23 perm. ; (3.28) 



K 



1 H ^ 



4 1 



12 



kik2k^ki 



i=l 



23 perm. . 



(3.29) 
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3.3 Summary of final results 

Here we summarize the final results from Sec. 13. H 13.21 and Appendix [Bl For the general 
single field inflation £(0, X), we define 



c: 



2 _ P,X 



Px + 2XPxx 
S = XP,x + 2X'P,xx , 

2 2 3 

X = X P^xx + -^X P^xxx 5 

/i = ^X^P,xx + 2X^P,xxx + lx^P,xxxx • (3.30) 

If any of /i/S, A^/S^, 1/cf > 1, the single field inflation generates a large primordial trispec- 
trum, whose leading terms are given by 

4 4 

(C') = (27r)9p|<53(^ k,,) J] ^ r(fci, ^2, fcs, ^4, fci2, fci4) , (3.31) 

j=l i=l i 

where T has the following six components: 

r^(A)V4(i-.)...(i-.)V..(|--V. 

'3A 1 /I 



+ (--7^ + 1 T,2+U-l)Te3. (3.32) 



The Tsi_s2,s3 are contributions from the scalar-exchange diagrams and are given in Ap- 
pendix [HI Tci^c2,c3 are contributions from the contact-interaction diagram and are given 
by fl3.27p - fl3.29p . For the most interesting cases, where any of yu/S, A^/S^, 1/c^ ^ 1, the 



first four terms in (13.320 are the leading contributions. So we have four shapes determined 
by three parameters jj A/S, 1/c^ and fi/T,. A large bispectrum necessarily implies a large 
trispectrum, because either 1/c^ or (A/S)^ is large. But the reverse is not necessarily true. 
One can in principle have a large /i/S but small 1/c^ and (A/E)^. 

To quantify the size (i.e., magnitude) of the non-Gaussianity for each shape, we define 
the following estimator t^i for each shape component, 

(C')component ^ {2nYP^5%y2 ^^'^79 , (3.33) 

limit ^ ' It 

i 

where the RT limit stands for the regular tetrahedron limit (fci = ^2 = ^3 = ^4 = ki2 = 
ki4 = k). The parameter ^atl is analogous to the /tvl parameter for bispectra. This definition 



^More generally, we have six shapes controlled by three parameters. However, the second line of p.32p are 
negligible unless /i/S, A^/E^, are all ^ 1, in which case the trispectra is only marginally large, ^ 0(1)- 
Also note that, for /i/S, A^/S^, ^ 1, the second line of (|3.32p docs not capture all the subleading 

contributions. 
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applies to both the cases of interest here, and the non-Gaussianities of the local form that we 
will discuss shortly. Unlike the convention in the bispectrum case where the normalization 
of Jnl is chosen according to the local form non-Gaussianity, here we conveniently choose 
the normalization of t^i according to fl3.33p . This is because, for the trispectra, even the 
local form has two different shapes. 

The size of non-Gaussianity for each shape in f l3.32p is then given by 

fj^L = 0.250 (^^^ , = 0.420^ (^^ - 1^ , = 0.305 (^^ - 1 

= 0.0352 - , = 0.0508 - 1 + 1^ , = 0.0503 (^^ - 1^ (3.34) 

For comparison, let us also look at the trispectrum of the local form. This is obtained 
from the ansatz in real space [39,40], 

C(x) = C. + pNL {C, - (O) + ^9NL {C! - 3(C,')C.) , (3.35) 

where (g is Gaussian and the shifts in the 2nd and 3rd terms are introduced to cancel 
the disconnected diagrams. Such a form constantly arises in multi-field models, where the 
large non-Gaussianities are converted from isocurvature modes at super-horizon scales. The 
resulting trispectrum is 

= InL^IocI + gNLTloc2 ■ (3.36) 

The two shapes are 

3 



'-loc2 



r,„, = ^(^ + 11 perm. I . (3.37) 

4 

100 

1=1 

where the 11 permutations includes fci3 A;i4 and 6 choices of picking two momenta such 
as ki and k2. The size of the trispectrum for each shape is 

t^l^ = 2.16/^^ , t%f = l.OSgML. (3.39) 

So again a large bispectrum implies a large trispectrum, but not reversely. 



4 Shapes of trispectra 

In this section, we investigate the shape of the trispectra. We take various limits of the 
shape functions T^i, Ts2, T^s and Td, and then compare among themselves, and with the 
local shapes Tiod and Tioc2- We will summarize the main results at the end of this section. 
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Figure 2: This figure illustrates the tetrahedron we consider. 



Before the discussion of the shape functions, we note that the arguments of the shape 
functions are six momenta ki, k2, k^, k^, ku, ki^. In order for these momenta to form a tetra- 
hedron (as in Fig. [2]), the following two conditions are required: 

Firstly, we define three angles at one vertex: 

cosfa) 



kj 


-f- 


■ k^ 




2fci/Ci4 






+ kl^^ - 


■ kl 




2/^2^14 




kf 


+ kl- 




2kik2 



cos(/5) 

These three angles should satisfy cos(a — 13) > cos(7) > cos(a + [3). This inequality is 
equivalent to 

1 - cos^(a) - cos^(/3) - cos^(7) + 2 cos(a) cos(/3) cos(7) > . (4.2) 

Secondly, the four momenta should satisfy all the triangle inequalities. We need 

ki + kA> ki4 , ki + k2> ki2 , k2 + k^> ku , 
kl + ku > ki , ki + ki2 > k2 , k2 + ku > k^ , 

k4 + ki4 > kl , k2 + ki2 > kl , ks + ku > k2 . (4.3) 
The last triangle inequality involving {k^, k^, ki2) is always satisfied given Eq. (14.21) and Eqs. 

dH). 
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We also would like to mention a symmetry in our trispectrum. As ki, k2, k^, are 
symmetric in our model, we have 

T{ki, k2, ks, k4, ki2, k^) = T{ki, k2, k^, k^, k^, k^) = T{ki, k^, k2, k^, k^, ku) , (4.4) 

and etc, where 

fci3 = |ki + ksl = ^kf + kl + kl + kl- kf^ - kf^ . (4.5) 

The first set of limits we would like to take is those involved in the consistency relations. 
There are two known consistency relations for the trispectra to satisfy. 

Firstly, we discuss the consistency relation in the squeezed limit. When one external 
momentum, say, k^ goes to zero, this mode can be treated as a classical background for the 
other modes, and the trispectrum should reduce to the product of a power spectrum and a 
running of bispectrum [26,37,38]: 

(CkiCk^CkaCk^) ~ -P(fc4)Tl^(CkiCk2Ck3) , (4.6) 

a ma 

where P(fc) = ^P^(/c), and P(^{k) is the dimensionless power spectrum. 

In our case, the leading order contribution to the trispectra scales as c"'^, or c~^A/S, or 
(A/S)^, or /i/S. However, RHS scales as c~^, or c^^A/S. In order that Eq. (14. 6 p holds, 
^4(CkiCk2Ck3Ck4) niust Vanish at the leading order in the ^4 — > limit. One can check that 
our result indeed vanish in this limit, Tsi 2,3,Tci —>■ 0{kl). 

Secondly, we check the folded limit, say ki2 —>■ 0. For the s-channel (in which the 
exchanged scalar carries the momentum ki2), the four-point function can be regarded as a 
pair of two-point functions modulated by the same classical background generated by the 
long wave mode ki2, and we have [28] 

(Ck.Ck.Ck^CkJ ~ {ns - l)'P(A:i)P(A;3)(C-k,,Ck,,) ■ (4.7) 

Note that the RHS takes the same shape as (I3.37p . Again, RHS scales as cj^ in our case. 
So in the ki2 limit, we expect fci2(CkiCk2Ck3Ck4) to vanish for the s-channel. For the 
t-, u- and the contact interaction channels, there are neither propagators that give rise to 
the pole behavior l/kf2 nor inverse Laplacians in our Lagrangian. Therefore in our case 
^i2(CkiCk2Ck3Ck4) trivially for these channels. 

One can check that our results indeed satisfy the condition. In fact, we have (CkiCk2Ck3Ck4) 
0{ki2) for the s-channel, so the pole behavior at ki2 = is cancelled more than enough to 
satisfy the condition. (Note that summing over all channels gives (C^) —>■ constant.) 

After checking the consistency relations, now we shall plot the shape functions. To do 
so, we shall take various limits to reduce the number of variables. We set the shape function 
to zero when the momenta do not form a tetrahedron. We consider the following cases: 
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1. Equilateral limit: ki = k2 = k^ = k^. In Fig. [21 we plot T^i, T52, T^^, T^, Tiod and Tioc2 
as functions of k^/ki and ku/ki. (We would like to remind the reader that unlike the 
first four shape functions, Ti^d and T/oc2 are not obtained in our model. We plot them 
for the purpose of comparison.) One observes that T/oci blows up at all boundaries. 
This feature can distinguish our shape functions from the local shape Tioci originated 
from the local fNL- 

2. Folded limit: ki2 = 0. (This limit is also related to the parallelogram limit, ki = ks, 
by the symmetry (14. 4p .) In this limit, ki = k2 and k-^ = k^. We plot T^i, Ts2, T53, Td 
and Tioc2 as functions of k^^/ki and ki^/ki in Fig. HI (Note that Tioci blows up in this 
limit). We assumed k^ < ki without losing generality. Note that Tioc2 does not vanish 
in the ^4 — > limit. This can be used to distinguish our shape functions from the local 
shape originated from g^i- 

3. Specialized planar limit: We take ki = k^ = kn^ and additionally the tetrahedron to 
be a planar quadrangle. In this limit, one can solve for ki2 from ( 14. 2^ : 



The minus sign solution can be related to another plus sign solution in the ki = k2 = ki^ 
limit through a symmetry discussed in Appendix O We will only consider the plus 
sign solution in our following discussion. We plot the shape functions as functions of 
k2/ki and k^/ki in Fig. [5l These figures illustrate two important distinctions between 
our shape functions and the local form shape functions. At the ^2 ^4 limit, we have 
^13 ~^ 0, so T/oci blows up, while the others are all finite. At the ^2 and ^4^0 
boundaries, our shapes functions vanish as 0{kl) and 0{kl) respectively, while Tiod 
and Tioc2 are non-vanishing. 

4. Near the double-squeezed limit: we consider the case where k^ = k4^ = ki2 and the 
tetrahedron is a planar quadrangle. We are interested in the behavior of the shape 
functions as k^ = k^i = ki2 — > 0, i.e. as the planar quadrangle is doubly squeezed. In 
this case, Eq. (14.21) takes the equal sign. One can solve for k2 from (14.21) . The solution 



is presented in Eq. (EI]). We plot T,i/(nti^*), T,2K]\k.i), T,.J{]\k.i), TdKUh), 



T;oci/(n ^i) ^"^^ Tioc2/iX\ ki) as functions of ^12/^1 and ku/ki in Fig. [61 To reduce the 
range of the plot, we only show the figures partially with k^ < ki. Note that in this fig- 
ure, we divided the shape functions by Y[ ki in order to have better distinction between 
contact-interaction and scalar-exchange contributions. Fig. [HI shows simultaneously 
the three differences among the four shapes T^i (~ ^^2,3), T^, Tiod and Tioc2- 1) In 
the double-squeezed limit, k^ = k^ ^ 0, the scalar-exchange contributions Tsi/{Ylki), 




(4.8) 
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Ts2/([l^i)y '^ss/iYl^i) nonzero and finite, and the contact-interaction Tci/{Y[ki) 
vanishes. As a comparison, the local form terms Tioci/{Y[ki) and Tioc2/{Y[^i) blow 
up. 2) In the folded limits, at the {k^/ki = l,ki4,/ki = 0) corner where fc^ — > 0, 
and close to the {k^/ki = l.ki^/ki = 2) area where ki^ 0, Tiod/ ([l^i) blows up. 
3) In the squeezed limit, at {k4/ki = l,ki4/ki = 1) where k2 ^ 0, the Tioci/iYl^i) 
and T;oc2/(n blow up. The last two behaviors have also appeared in the previous 
figures. 

In the second, third and fourth limits, the tetrahedron reduce to a planar quadrangle. 
We collectively denote this group of limits as the planar limit. This planar limit is of special 
importance, because one of the most important ways to probe trispectrum is the small 
(angular) scale CMB experiments. These experiments directly measure signals contributed 
mainly from the planar quadrangles. The more general plot for the planar limit is presented 
in Appendix O We can see that while very different from the two local shapes, the three 
shapes T^i, T52 and T^s are overall similar. Of course like in the bispectrum case, we can 
tune the parameters to subtract out the similarities and form new bases for the shapes. 

We end this section by emphasizing a couple of important points: 

• The equilateral trispectra forms: The scalar-exchange contributions T5i^2,3 and the 
contact-interaction contribution Td are similar at most regions, but can be distin- 
guished in the double-squeezed limit (e.g. ^3 = A;4 — > 0), where the two kinds of forms 
approach zero at different speeds, Tsi^2,3 ^^(^3); ^ci ^ ^^(^3)- Within the scalar- 
exchange contributions, the three shapes Tgi, Ts2, Tgs are very similar overall, having 
only small differences]^ 

For the purpose of data analyses, one can then use the following two representative 
forms for the "equilateral trispectra". One is Td, given in f l3.27p . This ansatz can be 
used to represent all four leading shapes at most regions. For a refined data analysis, 
for example to distinguish shapes in the double-squeezed limit, one can add another 
form Tsi, given in (IB. 31) and ( IB. 41) . This ansatz represents very well the three scalar- 



exchange contributions T5i^2,3- The first ansatz is factorizable (in terms of the six 
variables ^i, 2,3,4, ^12, ^14) by introducing an integral = (l/r(n)) Jg°°t"~^e~^* [4]; 

while the second ansatz cannot be easily factorized due to the presence of ^13 given by 

<KB. 



Distinguishing between the equilateral and local forms: In the following limits, the 
equilateral and local forms behave very differently. At the folded limit (e.g. ki2 



^For example, in Fig. [3] or 21 if we look at the double folded limit, ki2 and fci4 0, T32 and Ts3 go 
from positive to negative {Ts2 ^ —0.066 and Ts3 —0.030), while Tgi remains positive (T51 — + 0.092). 
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Figure 3: In this group of figures, we consider the equilateral limit ki = k2 = = k^, and 
plot Tsi, Ts2, Ts3, Tci, Tioci and Tioc2, respectively, as functions of ki2/ki and ku/ki. Note 
that Tioci blows up when ki2 <^ ki and ki^ <^ ki. Tiod also blows up in the other boundary, 
because this boundary corresponds to ki^ <^ ki. So Tioci is distinguishable from all other 
shapes in this limit. We also note that and Tioc2 are both independent of ki2 and k^. 
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1.0 



Figure 4: In this group of figures, we consider the folded hmit ki2 = 0, and plot T^i, Ts2, 
Ts3, Tci and Tioc2, respectively, as functions of k^/ki and k^/ki. Tiod blows up in this limit. 
Note that when — > 0, all shape functions except Tioci and Tioc2 vanish. 
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Figure 5: In this group of figures, we consider the speciahzed planar hmit with ki = = ku, 
and plot T51, Ts2, T53, Td, Tiod and Tioc2, respectively, as functions of ^2/^1 and k^/ki. Again, 
in the k2 ^ Q 01 k^^ ^ Q limit, our shape functions vanish as 0{k2) and 0{kl) respectively. 
This is different from that of the local shape. Tiod blows up when /c2 — > k^. This is because 
in this limit, kis —>■ 0. 
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Figure 6: In this group of figures, we look at the shapes near the double squeezed limit: we 
consider the case where ^3 = ^4 = ki2 and the tetrahedron is a planar quadrangle. We plot 

Tsi/ (Ilt=i ^i)' ^^2/ (n ^i)' Tss/ (n ^i)^ ^ci/(n ^t)^ Tiod/ (n ^i) and T10C2/ (n ^i)^ respectively, 
as functions of ki2/ki and ^14/^1. Note that, taking the double-squeezed limit ^4 —>■ 0, 
the scalar-exchange contributions Tsi/{Y[ki), Ts2/ ([l^i)^ ^ss/dl^*) nonzero and finite, 
and the contact-interaction Tci/{Ylki) vanishes. As a comparison, the local form terms 
Tioci/iYl^i) T10C2/ (Yl^i) blow up. The different behaviors in the folded and squeezed 
limit can also been seen from this figure (see the main text for details). 
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0), Tioci generically blows up, while the four equilateral shapes and Tioc2 approach 
constants. At the squeezed limit (e.g. ^4 — > 0), the four equilateral shapes all vanish 
as 0{kl), while the local forms Tiod and Tioc2 do not. 



5 Examples 



For DBI inflation [41-46], P = -/(0)-i/r^2X7(0) + /(0)-i - F(0), 

The dominant contribution come from the scalar-exchange terms %i,s2,s3 and one contact- 
interaction term T^, which are of order l/c^. The 7^2 and 7^3 are of order 1/c^, so belong 
to subleading contributions. Therefore the shape function defined in f l3.32p takes the form 

r^^'^{^L + I^ + T,^_Td)l., (5.2) 



4 2 c 



According to the definition (lOHD . t^f ^ 0.542/c^ 

For k-inflation [47-49], we look at the example P ~ {-X + X'^)/(j?, 

A 2X l-cl /i X l-cl 
< 1 — = = 1 tl = = i . (5.3) 

' ' S -1 + 6X 2 ' S -1 + 6X 4 ^ ' 

The only dominant term is one of the scalar-exchanging terms 7^3, the others all belong to 
subleading contributions. So 

, (5.4) 

and ^ 0.305/c^. 

As mentioned in the previous section, there are some differences among the shapes T^i, Ts2 
and Ts3, and especially between them and Td. These differences may be used to distinguish 
some special models within this class. But unfortunately, for the above two examples, after 
summing over all contributions for DBI inflation, the trispectrum difference between the 
DBI inflation and this specific k-inflation example becomes smaller, and we do not find any 
features that can very sharply distinguish them. This is because the four leading shapes 
Tsi and Td are similar in most regions and in the discriminating double-squeezed limit, the 
trispectra in both examples take the form Tgi (which are similar among themselves) since 
Tci vanishes faster. 
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6 Non-Bunch-Davies vacuum 



We now study the shape of trispectrum if the initial state of inflation deviates from the 
standard Bunch-Davies vacuum of de Sitter space. This is an interesting question because 
short distance physics may give rise to such deviation [52,53], and one might argue whether 
its effects on the power spectrum of the CMB are observable [54 ^. The effects of the non- 
Bunch-Davies vacuum on the bispectrum have been studied in [9,11,12], where it was found 
that the non-Gaussianities are boosted in the folded triangle limit (e.g. ki + k2 — k-^ ^ 0) . 
A general vacuum state for the fluctuation of the inflaton during inflation can be written 

as 

Uk = u{k, t) = -^=£==(C+(1 + zA;c,r)e-*^^=^ + C_(l - tkcsT)e"'''n ■ (6-1) 

Here a small and non-zero C_ parametrizes a deviation from the standard Bunch-Davies 
vacuum which has (7+ = 1, C_ = 0. We consider the corrections of a non-zero C_ to the 
leading shapes assuming C_ is small, we keep only terms up to linear order in C_. Similar 
to the case of the bispectrum [9], the first sub- leading corrections come from replacing one 
of the M(r, k)'s with their C_ components, and since the correction from M(0,k) only has 
the same shape as that of the Bunch-Davies vacuum, we only need to consider the case that 
u{t, k) comes from the interacting Hamiltonian, where r is not zero. In the following we 
discuss the contact-interaction diagram and the scalar-exchange diagram respectively. 

For the contact-interaction diagram, the corrections consist of four terms from replacing 
ki with —ki in the shape for Bunch-Davies vacuum. For the leading shape Td we denote the 
correction as and we find 



f,i = 36 Re(C_) JJA;,^ 



=1 



1 1 

+ 



(/ci + k2 + ks- kif {ki + k2-kz + k^) 



1 1 
+ ; ; + 



{ki - k2 + k^ + kif i-ki + k2 + ks + ki)\ 



(6.2) 



Then let us consider scalar-exchange diagram. For illustration we only consider the Tgi 
term. The calculations for Ts2 and Tss are similar but more complicated. Now we have six 
Uk modes from the interaction Hamiltonian. Replacing each mode with its C_ component 
gives rise to six terms in the corrections. They correspond to replacing ki with —ki, or one 
of the two /ci2's with —ki2 in the calculations for the Bunch-Davies vacuum. The corrections 



^The choice of initial state is often discussed in the context of trans-Planckian effects [55] though the 
issue has more general applicability. See e.g. [56] for a review and references, and [57] for a discussion of 
how to capture the initial state effects in terms of a boundary effective field theory. 



21 



to Eq. (1R3|) are 

9 Re(C_ 



k-^ }xi<2 /c^ ^12 



E 



1 



(/Cl + k2- ki2f {h + /C4 + fcl2)^ 

+23 perm. 



ih + h + ki2f {ki + k2 + ku) 



+ 



(/ci + k2 + k^2f ik3 + ki- ki2f 



(6.3) 



In Eq. ( ]B.4[) . the two ki2 cancelled in the final expression for the Bunch-Davies vacuum, 
so we have to recover them in the calculations. Denoting M = k^, + k^ + ki2 and K = 
^1 + ^2 + ^3 + ^4, we find the corrections 



9 Re{CJ] 



4 



k-^ ^12 



+ 



(ks- 
1 



- k,2)' 



6(A;3 + ki- ki2) 



+ 



{K - 2ku) 
3M 



1 

+ 



6M2 3M 1 



3{k3 + k^ — k 



12) 



+ 



1 



2k 



+ 



1 



12J 



{K - 2k,2] 
+ 23 perm. 



(6.4) 



XK + 2k^2f {K + 2k^2Y {K + 2k^2f. 
To summarize, the correction T51 to Tgi is the summation of Eqs. (16.31) and (16. 4p . 

We can look for the analogue of the folded triangle limit discovered in the study of 
bispectrum where the corrections due to deviation from the Bunch-Davies vacuum diverge. 
Here we see when any one triangle, e.g., (fci, k2-, ^12), in the momentum tetrahedron becomes 
folded, the corrections to Tsi (see (16. 3 p and (16. 4p ) become divergent. Furthermore, when ki + 
^2 + ^3 — ^4 = 0, the two triangles (/ci, ^2, /C12) and (^3, ^4, ^12) become folded simultaneously, 
and the correction (16.21) to Td also diverges. We will refer to such configurations as the 
folded sub-triangle configurations,^ 

As discussed in [9], these divergences are artificial, and do not correspond to real infinities 
in observables. Rather, the divergences appear because it is not realistic to assume a non- 
standard vacuum to exist in the infinite past. A cutoff on momenta should be imposed at 
the same time when a non-Bunch-Davies vacuum is considered. 

We would like to point out two interesting aspects of the effects of the non-Bunch-Davies 
vacuum on trispectra, and more generally on higher point functions. 

Let us first look at the regions away from the folded sub-triangle configurations. In the 
regular tetrahedron limit, in terms of (I3.33p . the corresponding t^L for the non-Bunch-Davies 
contribution are 

,2\ /x\2 



t 



cl 
NL 



4.50 Re(C_) ( I - ^ 



t 



si 
NL 



401 Re(C_) ( ^ 



(6.5) 



^Duc to permutations, the three momenta do not have to be next to each other in terms of Fig. [21 for 
example, fci, fca, fcia. 
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Figure 7: In this group of figures, we plot the Tci/Re(C_) (the left column) and Tsi/Re{C^) 
(the right column) in the equilateral limit, specialized planar limit, and near double squeezed 
limit (we plot Tci/[Re{C-)Iliki] and Tsi/[Re{C-)Iliki] in near double squeezed limit) respec- 
tively. Note that, in order to show clearly the locations of the divergence, in some figures 
we have taken the cutoffs of the z-axes to be extremely large. 
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Note that t^j;^/Re(C_) is 128 times larger than t^^; t^j;^/Re(C_) is about 1600 times larger 
than t^^. These large numbers arise because some plus signs become minus signs in the 
denominators, and there are also more terms to consider in the non-Bunch-Davies case. So 
in the context of general single field infiation, even if Re(C_) is as small as one part in 
one thousand, it becomes important phenomenologically. Generalize this to higher point 
functions, we see that, no matter how small the C_ is, there always exists a high point 
function beyond which the contribution from the C_ component becomes comparable to 
that from the C+ component. For such functions, we should use the whole wavefunction 
(16. ip to get the correct shapes instead of treating the C_ component as a correction, even 
away from the folded sub-triangle limit. Therefore generally speaking, higher point function 
is a more sensitive probe to the non-Bunch-Davies component. However, on the other hand, 
higher point functions contribute less to the total non-Gaussianities and will be more difficult 
to measure experimentally. 

We next look at the region near the folded sub-triangle limits. Because the denominators 
here have larger powers than those in 3pt, the corrections grow faster as we approach the 
folded sub-triangle limit. This also indicates that the trispectrum, and more generally higher 
point functions, is a nice probe of the non-Bunch-Davies vacuum. But, on the other hand, 
because the higher point function has a larger momentum phase space (three more dimensions 
here in trispectra) than the 3pt, the phase space for the folded limits becomes relatively 
smaller. 

It will be interesting to see how the two factors in each of the above two aspects play out 
in the data analyses. 

In Fig. [71 we plot Tci/Re(C_) and Tsi/Re(C_) in the equilateral limit, specialized planar 
limit, and near double squeezed limit {Tci/[Re{C-)Iliki] and Tsi/[Re{C-)Iliki] in this case) 
respectively. 

7 Conclusion 

To conclude, we have calculated the leading order trispectra for general single field infiation. 
As in the case of bispectra, the trispectra turns out to be of "equilateral shape" in gen- 
eral single field infiation. Compared with the local shape trispectra, the equilateral shape 
trispectra has not been extensively investigated in the literature. It is clear that there are a 
lot of work worthy to be done on this topic in the future. Directions for future work include: 

• It is useful to extend our calculation to more general cases. We have focused here on 
the leading order contribution and single field infiation. It is interesting to generalize 
this calculation to next-to-leading order which may be also potentially observable, or 
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multifield inflation [58]. It is also useful to perform a unified analysis for general single 
field inflation and slow roll inflation. 

• The shape of the trispectra is much more complicated compared with that of the 
bispectra. In our paper, we have obtained a lot of features of the shape functions by 
taking various limits. However, it is still a challenge to find improved representations 
to understand the shape functions. For example, one can find new bases for the shapes 
by tuning parameters to subtract out the similarities. Also, we focus on the planar 
limit in plotting the figures (except for Fig. [3l), because the planar limit is of special 
importance for the CMB data analysis. It is interesting to investigate the non-planar 
parameter region in more details for the large scale structure and 21-cm line surveys. 

• In the discussion of non-Bunch-Davies vacuum, we have calculated contributions from 
the two representative shapes. However, as we have seen in Section 6, trispectra is a 
powerful probe for non-Bunch-Davies vacuum. Even one part in 10^ deviation from 
the Bunch-Davies vacuum could lead to an order one correction to the trispectra. So 
it is valuable to perform the full calculation for the non-Bunch-Davies vacuum, and to 
study the effects of the cutoff. 

• Most importantly, one would like to apply these shape functions to data analyses and 
see how they are constrained. 

Note added: On the day this work appeared on the arXiv, the paper [59] was also submitted 
to the arXiv, which overlaps with our Sec. 3 and 5. 
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A Commutator form of the in-in formalism 

In Sec. 13.11 we have used the original definition of the in-in formalism (13. Sp and (13. 9p to 
calculate the correlation function due to the scalar-exchange diagram. Another equivalent 
and commonly-used form is written in terms of the nested commutators. For the diagrams 



that we considered in Sec. 13. H it takes the following form, 



(Od-/ dt' [ dt"{0\[[Cj{t),Hj{t')],Hj{t")]\0) . (A.l) 

J to J to 



The main difference is that now the first term in (13. 9p is separated into two integrals, each 
has a time- or anti-time-ordered double integration. 



'dt' f dt" (0| Hj{t") Cf{t) Hj{t') |0) + f dt' f dt" (0| Hj{t') Cf{t) Hj{t") |0) . (A.2) 

'to J to -J to J to 

If one uses this form, besides the fact that the algebra becomes much more complicated due 
to the time-ordered double integrations, one also encounters spurious divergences at special 
momentum configurations, such as ki + k2 — — = 0, for each of the two terms in ( lA.2p . 
These divergences can be seen, after some complicated algebra, to cancel each other once 
the two terms are summed up [28,50,51]. Therefore using the form of the first term in 
(13. 9p is both algebraically simpler and free of spurious divergences. This conclusion can be 
generalized to the more nested terms. 

For example, using formula ( lA.ip we can get 

(CV - 1^2 ^, {3{KU + Kl,) + {l2kl,-9K,^,K,,,) 

» ^ ^(_)^(+) (2A;i2 - i^(-) + K^+)) 

X (^(-) + ^(+)) + 8^(-)^(+) i^U + ^(+)) - 6A:i2 [2 (Kf_) - irf^)) 
-3K^-)K^+) - } + 23 perm. , (A.3) 

with = fci + ^2 — /cs — A;4 and = /ci + A;2 + fc3 + fc4. Although the terms before 23 perm, 
are not equivalent to the one in ( 1B.3P plus ( lB.4p . however, after including the permutation 



terms and performing lots of lengthy but straightforward calculations, one can find that the 
two expressions are the same. 
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B Details on the scalar-exchange diagram 

In this Appendix, we give the details of the scalar-exchange diagram. We denote 

ki2 = ki + k2, M = k3 + h + ki2, K = ki + k2 + h + . (B.l) 
The following are the various contributions to the trispectrum form factor T defined as 

4 4 

(C") = {27,f Pl5\Y, k.) n p ^(^1' ^3, ku, ku) . (B.2) 

i=i i=i « 

The interaction Hamiltonian has two components, fl3.6l) and (13.71) . So there are four different 
combinations. 



B.l The component {('^) 



4\ 

aa 



This component is given in (13.151) and (I3.17p . The contribution from the first term of 
Eq. ([33]): 

The contribution from the second and third terms of Eq. (13. 9p : 
B.2 The component {C^)ab 

There are two sub-diagrams contributing to this case. In the first sub-diagram, the exchanged 
scalar propagator is due to the contraction between two C's; in the second, between ( and 
(. The result of these two sub-diagrams for the first term of Eq. (13.91) . and the second and 
third term of Eq. (13. 9p follows. 

The contribution from the first term: 

• The first sub- diagram: 

• The second sub-diagram: 

-^^(-2-^) (ki2 ■ k4) ^'f^%, ^, F{ku, k„ M) + 23 perm. . (B.6) 

16S Vc^ J ki2 [ki + k2 + ku)'^ 
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The contribution from the second and third terms: 

• The first sub- diagram: 

- (^^ - (k3 ■ k^)klklh2 Gabih, h) + 23 perm. . (B.7) 

• The second sub-diagram: 
-^^ f ^-ij (ki2-k4)^^ Gab{ki2,k,) + 23peTm. . (B.8) 
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B.3 The component {C^)ba 

Similar to the (C^)afe case, here we also have two terms, each term includes two sub-diagrams: 

The contribution from the first term: 
Same as ([RS]) and flRGl) . 



The contribution from the second and third terms: 
• The first sub- diagram: 



- (^^ - (ki ■ k,)klklk^2 GUki, k2) + 23 perm. . (B.9) 
The second sub-diagram: 

^ ^ / 1 ) (k2 ■ ki2)^r^ GU-ki2, k2) + 23 perm. . (B.IO) 

12 



8 s V c? r ^ k 



B.4 The component {C,^)bb 

In this case, we have four sub-diagrams for each term. In the first sub-diagram, the exchanged 
scalar propagator is due to the contraction between two C's; in the second, between ^ and 
in the third, between C, and in the fourth, between two C's. The result of these four 
sub-diagrams for the first term of Eq. (13. 9p . and the second and third term of Eq. (13.91) are 
as follows. 

The contribution from the first term: 
• The first sub-diagram: 

2 

The second and third sub-diagrams: 



\ (ki ■ k2)(k3 ■ k4)A;i2 F(A;i, ^2, ki + k2 + ku)F{k3, k^, M) + 23 perm. . (B.ll) 



(ki ■ k2)(ki2 ■ k4)|^ F(A;i, fc2, fci + ^2 + ku)F{ki2, k,, M) + 23 perm. (B.12) 
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• The fourth sub-diagram: 



~ ^ ~ ■ ^^^^^'^ ' ^^^^ ^1 + ^2 + h2)F{h2, h, M) + 23 pern<.B.13) 

The contribution from the second and third terms: 
• The first sub-diagram: 

26 

The second sub-diagram: 

(ki ■ k2)(ki2 ■ k4)^ 6-^6(^1, ^2, A;i2, k^) + 23 perm. . (B.15) 



\ (ki • k2)(k3 • k4)A;i2 A;2, A;3, A;4) + 23 perm. . (B.14) 



C2 



■'s 



12 



• The third sub-diagram: 



^f^-lj (ki2-k2)(k3-k4)^ G';,6(-A;i2,A;2,A;3,A;4) + 23perm. . (B.16) 
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• The fourth sub-diagram: 

~ ^ " ' ' '^"''^ <^;'''(-^i2, k2, ki2, k^) + 23 perm. . (B.17) 

The function F, Gab, G^a, G^b are defined as follows: 

F{ai,a2,m) 

= — [2q;iq;2 + {ai + a2)m + rri^] , (B.18) 

Gab{ai, a2) 
= [2aiOC2 + (ai + a2)M + M'] 

3 12 
+ ^^^^^'^ + ("1 + "2)M] + ;^^^"i"2 , (B.19) 

^ Wk^ M^^"^ + «2 + M) + ^ [2aia2 + 2{a, + a2)M + M^] 



^2, "3, "4) 
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+ 



+ 



+ 



1 
2 
6 



[2a3a4 + (as + a4)M + M^] 

4 

2a3a4(Q;i + 0:2) + (2a;3a;4 + (ai + 02) ("3 + «4)) M + ^ a^M 
4 

2 JJ^ «i + (2a3a4(ai + 0:2) + aiQ;2(a;3 + 04)) M + ajajM^ 

i=l i<j 
i=l / \ i=l / i=l 



(B.21) 



Note that in Gab-, Gba and Gf,b, the K and M are defined as K = fci + ^2 + ^3 + ^4 and 
M = A;3 + A;4 + ku, but not in terms of a^'s. 

To summarize we denote the overall contribution from the scalar-exchange diagrams as 

2 \ / 1 \ /I \ 2 



fB.22) 



where T^i is given by fIB.SP and flB.4p . T52 is given by flB.5l) - flB.10l) . is given by fIB.lip 

(EITD. 



C The planar limit of the trispectra 

In Sec. m we discussed various properties of the shape functions. As we have stated, the 
planar limit has special importance for CMB experiments. So in this appendix, we investigate 
the trispectra in more detail in the planar limit and perform a survey of parameters for the 
shape functions. 

In the planar limit, Eq. (14. 2 p takes the equal sign. One can solve ^2 from Eq. (14. 2p . 

. (_t.2 I 1.2 , t.2\ I L2 L2 I p L2 , i,2 l2 , 1.2 1.2 _ t.2 t.2 _ h.4 i l2 l2 
_ V ^1 I ^12 ^ ^3 ^ ^4J ^^l^s2 ^ ^12^14 ^ ^12^4 ^ ^^14^4 ^14^3 ^4 ^ ^3^4 j^q 



y2A;4 



where ksi and fcs2 are defined as 



/c^i = 2^{kiki + ki ■ k4)(A;iA;4 - ki ■ k, 
/c^2 = 2 a/ (A;3/C4 + k3 ■ k4)(fc3/c4 - k3 ■ k4 



4j , 



(C.2) 



We first take a close look at the ± sign in Eq. (IC.ip . The — sign and the + sign correspond 
to two different quadrangles. In Fig. [8], the — and + solutions correspond to the black (with 
edge ki) and blue (with edge qi) quadrangles respectively. The former has all internal angles 
< TT, while the latter has one > vr. The blue quadrangle can be transformed into another 
quadrangle belonging to the same class as the black one by the symmetry discussed in Eq. 
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Figure 8: The quadrangles in black and blue represents the solution fIC.ip with minus and 
plus sign respectively. The term with plus sign corresponds to the blue quadrangle, and can 
be transformed into another quadrangle corresponding to the minus solution by a symmetry 
discussed in Eq. f l4.4p . 

(14.41) . So without losing generality, we will only consider the — solution in the following 
discussion. 

To scan the parameter space. We plot T,!, Ts2, Tgs, Td, Tioci and Tioc2 as functions of 
ku/ki and ki^/ ki for different values of ks/ki and k^/ ki in FigjOl The momenta {ks/ki, k^/ ki) 
take the following values in the 4 figures in each group. 

{k^/ki, kjk,) = {(0.6, 0.6), (0.6, 1.0), (1.0, 0.6), (1.0, 1.0)} . (C.3) 

We assume ki > /i;2, ^s, ^4 in the plot without losing generality. Note that when k^ = k^ = 
ki2 = ^14 = ki (the center of the last figure in each group), T^i, Ts2, T^^ and vanishes 
because in this case k2 = 0. We can see from these graphs that the shapes of T51, Ts2 and 
Ts3 are overall very similar. 
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Figure 9: In the six rows, we plot T^i, Ts2, T'ss, T'ci, Tiod and T/oc2 respectively as functions 
of ku/ki and ku/ki in the planar limit. Within each row, the momenta configuration is 
{k3/ki,kjki) = {(0.6,0.6),(0.6, 1.0),(1.0,0.6),(1.0, 1.0)} respectively in the four columns, 
as given in (IC.Sp . 
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